function output_dat(A,name,n)
f=fopen(strcat(name,'_',num2str(n),'.dat'),'w');
for i=1:size(A,1)
    fprintf(f,'%e \n',A(i,:));
end
fclose(f);
end